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ABSTRACT 

We have observed the central region of the nearby starburst galaxy NGC 253 
with the Kyoto Tridimensional Spectrograph II (Kyoto3DII) Fabry-Perot mode 
in order to investigate the properties of its galactic wind. Since this galaxy has 
a large inclination, it is easy to observe its galactic wind. We produced the Ha, 
[N II]A6583, and [S II]AA6716,6731 images, as well as those line ratio maps. The 
[N II] /Ha ratio in the galactic wind region is larger than those in H II regions in 
the galactic disk. The [N II] /Ha ratio in the southeastern filament, a part of the 
galactic wind, is the largest and reaches about 1.5. These large [N II] /Ha ratios 
are explained by shock ionization/excitation. Using the [S II] /Ha ratio map, 
we spatially separate the galactic wind region from the starburst region. The 
kinetic energy of the galactic wind can be sufficiently supplied by supernovae in 
a starburst region in the galactic center. The shape of the galactic wind and the 
line ratio maps are non-axisymmetric about the galactic minor axis, which is also 
seen in M82. In the [N II]A6583/[S II]AA6716,6731 map, the positions with large 
ratios coincide with the positions of star clusters found in the Hubble Space Tele- 
scope (HST) observation. This means that intense star formation causes strong 
nitrogen enrichment in these regions. Our unique data of the line ratio maps 
including [S II] lines have demonstrated their effectiveness for clearly distinguish- 
ing between shocked gas regions and starburst regions, determining the extent 
of galactic wind and its mass and kinetic energy, and discovering regions with 
enhanced nitrogen abundance. 
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Subject headings: galaxies: individual(NGC 253) — galaxies: starburst — ISM: 
abundances — ISM: jets and outflows 



INTRODUCTION 



Galactic scale outflows (galactic winds) from galact ic disks are common phenomena in 



the u niverse: they are found both in nearby galaxies (e.g., Strickl and et a 



20061 ) and in high-z galaxies (e.g., IShapley et al.l 120031 ; iTapken et al 



2004l : lTullmann et al. 
Their main 



2003). 



sources of energy are stellar winds, supernovae, or active galactic nuclei (AGNs). Their 
asy mmetric and comp lex structures (e.g., M82; IShopbell fc Bland-Hawthornl Il998l . NGC 



253; iSugai et al.l 120031 ) are investigated by numerical simulations (e.g., iCooper et al 



Galac tic winds are considered to impact on galaxies and their environments. iKobayashi et al. 



2008). 



( 120071 ) takes galaxy-scale outflows into their semi-analytical model of galaxy formation 
as a parameter, which r eproduces t he ob served Lya luminosity functions of Lya emit- 
ters in high-z. Recently, ISato et al.l (120081 ) has suggested from their multi-wavelength ob- 
servation data that galactic winds may carry away interstellar gas in host galaxies and 
may quench star formation. Galactic winds from dwarf galaxies with gas masses of a 
few xlO 8 M are suggested to be the m ajor pollutants of the in t ergalactic med i um (I GM) 
through an analytical / n umeri cal model (IFerrara &: Tolstoy! l2000l ). IVeilleux et al.l (120051 ) and 
Bland-Hawthorn et al.l (120071 ). in reviews of the galactic winds, summarize the properties 
and significance of galactic winds. 

Despite the importance of galactic winds, there are few galaxies where the basic physical 
values (e.g., mass and kinetic energy) of the galactic winds were derived from spatially 
resolved data. It is important to analyse spatially resolved data because galactic winds have 
extremely complex st ructures: e.g., non-axisymmetry a bout the minor axis of a galacti c 
disk is found in M82 Jshopbell fc Bland-Hawtho7nlll998h and NGC 253 JSueai et al.l l2003h . 
The Fabry- Perot Interferometer, which enables us to subtract the continuum precisely and to 
produce pure line images, is a powerful meth od to observe gala ctic winds. The galactic winds 
in M82 (IShopbell fc Bland-HawthoTr] Il998h and NGC 1482 JVeilleux fc Rupke! l2002h were 
observed with Fabry-Perot instruments, and their masses, kinetic energies, and excitation 
sources of them were estimated from the [N II]A6583/Ha line ratio maps. In spite of the fact 
that this method is especially effective, there were only three examples where the properties 
of galactic winds in edge-on galaxies were discussed in detail from line ratio maps. 

NGC 253 is one of the most famous g alaxies beca use it is a nearby (D=3.5 Mpc: 



Rekola et al.ll2005l ) edge-on (inclination=78°: |Pencelll980f ) starburst galaxy. Details of thi s 



galaxy have been investigated in many wavelengths such as X-ray (IStrickland et al.l I200CH ) 



- 3- 



optic al (ISchulz fe Wegnerlll992l; Watson et al.lll996[). infrared (ISams et al.lll9 94; Su gai et al. 



1997). In NGC 253, the galactic wind was found in Ha (e.g. . McCarthy et al.lll987l ; ISchulz fe Wegner 



2003 ). sub-mm ( Sakamoto et al. 2006 ). and radio ( Turner fe Hql985 ; Ulvestad fe Antonucci 



19921 ) and in X-ray (e.g. JStrickland et al.ll2000l ). and the properties of the galactic wind were 



discussed in th em. The expansion vel o city is estimated t hroug h slit spectroscopic observa- 
tion in optical (ISchulz fe Wegnerlll992l ). IStrickland et al.l (120001 ) discuss the spatial structure 
of the outflow from the azimuthal surface brightness profiles of the X-ray and the Ha. 

We have observed the central region of NGC 25 3 with the Fabry- Perot (FP) m ode of the 
Kyoto Tridimensional Spectrograph II (Kyoto3DII: ISugai et al.ll2000l . 120021 . 120041 ) mounted 
on the Subaru Telescope. The Kyoto3DII is a unique instrument that enables us to carry 
out the optical Fabry-Perot observations on an 8-meter class telescope and to observe weak 
emission lines. The observations were carried out during the test observations of the Ky- 
oto3DII, and this is the first result of the Kyoto3DII FP mode. From these data, we have 
obtained the [N II]A6583/Ha, [S II]AA6716,6731/Ha, and [N II] /[S II] line ratio maps. This 
unique set of line ratio maps allows us to clearly distinguish between shocked gas regions and 
starburst regions, to determine the extent of galactic wind as well as its mass and kinetic 
energy, and to discover regions with enhanced nitrogen abundance. 

We describe our observations and data reduction in Section [2j We show the Ha and 
continuum images and the line ratio maps, and compare our images with the X-ray and the 
radio images in Section [3J We discuss the properties of the galactic wind and stellar clusters 
from the continuum and Ha images and the line ratio maps in Section |H We summarize our 
conclusions in Section [5j 



OBSERVATIONS AND DATA REDUCTION 



We have o bserved the central region of N GC 253 with the Fabry-Perot (FP) mode of 
the Kyoto3DII JSugai et al.l boool . I2OO2I I2004T ) on 28 August 2002 during the test observa- 
tion of the Kyoto3DII. It was mounted on the Cassegrain focus of the Subaru Telescope 
at Mauna Kea. This mode uses an ET-50 etalon made by Queensgate Instruments. The 
field of view is ~ 1'.9 x 1'.9, with the spatial sampling of 0".112 after 2x2 on-chip bin- 
ning. The spatial resolution was 1".0 and the spectral resolution was 21 A. For each of the 
Ha+[N II]AA6548,6583 and [S II]AA6716,6731 sets, we obtained 5 on-band frames and 1 off- 
band frame. The central wavelengths of the transmission curve of etalon for the Ha+[N II] 
set were 6556 A, 6568 A, 6580 A, 6592 A, 6604 A(on-band), and 6502 A(off-band), while those 
of the [S II]AA6716,6731 set were 6717 A, 6729 A, 6741 A, 6754 A, 6766 A(on-band), and 
6650 A(off-band). The exposure time of each frame was 240 seconds. The bias subtraction 
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and the flat fielding have been performed for the target object frames and the standard star 
frames, and then the flux calibration has been done. For reducing the observing time, flat 
and standard star frames were not obtained at all the observed wavelengths. Therefore, we 
have used the flat frames of the nearest wavelength for the flat fielding and the interpolated 
data for the flux calibration. The flat flames in 6502 A, 6568 A, and 6650 A were similar to 
each other: ~ 20 % brighter at the frame center compared with at the edge. The flat flames 
in longer wavelengths, including in 6729 A, 6742 A, 6749 A, and 6757 A, had a different flat 
pattern: darker at the frame center. In the longest wavelength, it was ~ 30 % darker at the 
center. The resultant uncertainty caused by the flat fielding, including the interpolation, is 
negligible in shorter wavelengths, and within a few percent in longer wavelengths. We have 
removed cosmic rays, matched the spatial resolution to the worst frame, and subtracted the 
off-band frames from the on-band frames. Eventually, the frames including only emission 
lines were obtained. The sky subtraction has not been performed in each frame because 
objects are observed in the whole field of view. This affects the continuum image, but does 
not affect the line images because the sky emission has been subtracted through the process 
of the off-band subtraction. 

Because of the low spectral resolution, the Ha and the [N II]AA6548,6583 in the Ha+[N II] 
set, as well as the [S II]AA6716,6731 in the [S II] set, were blended. In order to decompose 
them, we have fitted these lines pixel by p ixel with Airy functio ns, which are the transmission 



curve of the Fabry- Perot Interferometer (IBland &: Tullylll989l ). Figure [T] shows examples of 



the fitting. The residual ratios were calculated by dividing the total residuals by the total 
flux of the Ha plus [Nil] AA6548, 6583, and they were 5% and 30% in good and poor exam- 
ples, respectively. In poor fitting, line decomposing may be inaccurate, and the flux ratio 
may have a large error. However, the fraction of poor examples was smaller than 1/100, and 
much smaller in the regions we discuss below, and we neglected them. The average residual 
ratio for a single pixel in the discussed regions was 15%. As free parameters, we took the 
fluxes of the Ha and [N II]A6583 in the Ha+[N II] set, and those of the [S II]A6716 and 
[S II]A6731 in the [S II] set. We assumed that the flux of the [N II]A6548 was one third of 
that of the [N II]A6583. For the fitting, we fixed the vel ocity center so that it corresponded 



to the heliocentric velocity of NGC 253 of 251 km s _1 (IWatson et al.lll996l ). We also fixed 
the intrinsic velocity dispersion as zero because it was much smaller than the instrumental 
profile. The main source of the uncertainties for the line ratio maps is not flux calibration 
but line fitting, because the sky transmittance was stable during the observation for NGC 
253. In order to estimate the uncertainties of the fitting, we considered pseudo-observation 
data: they were obtained with the same "observing" parameters as actual observations from 
the spectra which include maxi mum velocity variation, velocity dispersion, and line-splitting 
observed in slit-spectroscopy (jSchulz fe Wegnerl Il992l ). The uncertainties were estimated 
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from the difference between the ratio derived from the fitting and the input intrinsic ratio. 
The uncertainties of the [N II] /Ha and [S II] /Ha maps were smaller than 30% and 10%, 
respectively. The uncertainty of the absolute flux calibration was less than 20%, and the un- 
certainty of the Ha flux, including the line fitting, was 30%. As for the [S II]A6716/[S Il]A6731 
ratio, it can be largely affected by the velocity field. Therefore, the [S II] ratio map is not 
presented. There was a distortion aberration, as expected, in the frames and we have cor- 
rected it. The background noise level was 1.9 x 10~ 16 erg cm -2 s _1 arcsec -2 in the Ha set 
and 6.3 x 10 -17 erg cm -2 s -1 arcsec -2 in the [S II] set, respectively. We determined the ab- 
solute coordinate by comparing our FP continuum fr ame with the contin uum frame of HST 
Wide-Field Planetary Camera 2 (WFPC2) fF675W: Iwatson et al.l fl99~6l ) . The uncertainty 
of the coordinate between the FP and the H ST frames was e stimated to be 0".08 {la). We 
adopt a distance of 3.5 Mpc to NGC 253 jRekola et alil2005h . 



RESULTS 



3.1. The continuum and Ha images 



Figure [2] (a) shows the continuum image of central region of NGC 253. Figures [3] (a), [2] 
(b), and [2] (c) show the whole field of view, the central region (deep), and the central region 
(shallow) Ha images, respectively. The continuum image is obtained from the average flux 
of the two off-band frames. In the central region, two bright regions are found in the Ha 
and continuum emission. One is found at 4"southwest from the center of NGC 253, while 
the other is found at 5"southeast from the center. Hereafter, we call the former starburst 
region A, and the latter starburst region B. A filamentary structure is seen to the northeast 
of the center and we call it filament C. In regions farther from the center, there are two 
bright regions in the Ha image. One is found at ~30"east from the center, and the other is 
found at ~60"southwest from the center. We call the former H II region A, and the latter 
H II region B. We analyse the characteristics of these regions in Section 13.21 

In the Ha and [N II] images, two faint filaments are seen to the southeast of the center. 
They were found in the narrowband Ha+[N II] image ([McCarthy et al.lll987l : the Ha and 
[N II] lines are blended in their image). These filamentary structures are obse rved at the edges 



of th e outflow when we observe the conical shell-like optically-thin gas (e.g.. lStrickland et al. 



2 OOOl ; ISugai et al.ll2003l ). The positions of these fil amentary structur es coincide with those in 
the near-infrared H2 emission (features A and B: ISugai et al.l 120031 ). Hereafter, we call the 
southwest one filament A and the southeast one filament B. We refer to these filamentary 
structures as "filaments", although we consider that they are not filaments but actually are 
a part of conical outflow. The north-side features were also found in the H2 emission, but 
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not in the Ha probably due to stron g extinction. The posi tions of filamentary structures 
seen in the Ha and in the soft X-ray ({Strickland et al.ll2000l ) in the galactic wind are close 
to each other (Figures [2] (h) and [3] (b)). This indicates that the relatively warm (~ 10 4 K) 
gas exists closely to the hot (~ 10 6 K) gas. A closer inspection, however, reveals that the 
X-ray filaments tend to exist where the Ha emissions are locally weak and tend to exist 
in inner regions compared with the Ha filaments: the separation between the Ha and the 
X-ray filaments is ~ 70 pc. This represents that the inner part of the galactic wind has a 
higher temperature than the outer part of it. This relation between the spatial distributions 
of t he Ha and the X-ray h as been found through our accurate positioning of these two maps 
(cf. IStrickland et al.ll2000i ). The spatial resolution of our Ha image (=1".0) also is slightly 
better than that obtained by Strickland et al. 200ol (=1".5). 



Starburst region B is as bright as starburst region A in the Ha a nd the cont i nuum 
at ~ 6600 A, while t he counterpart is not found in the near- infrared (ISams et al.l Il994l ) 
or radio observations (jUlvestad fc Antonuccil 119971 ). This region is diffuse in the optical 
broadband images r etrieved from HST arch ives. A comparison of our Ha map and long-slit 
spectroscopic data (jSchulz fc Wegnerl Il992l ) has revealed that starburst region B has the 
sharp local-minimum velocity of -100 km s _1 . 



3.2. Line ratio maps 



Figures[2](d) and 121(e) show the [N II]A6583/Ha and [S II]AA6716,6731/Ha ratio maps, 
respectively. The [N II] /Ha line ratio varies widely among regions. The [N II] /Ha ratios of 
filaments A and B are larger than the other regions. The ratio is particularly the largest at 
the tip of filament B and reaches about 1.5. As leaving from the galactic center along the 
minor axis, the Ha flux of the galactic wind becomes weaker while its [N II] /Ha ratio tends 
to become larger. The [N II] /Ha ratios of starburst regions A and B are ~ 0.5, which are 
local minima in the central region and the galactic wind. The [N II] /Ha ratios of H II regions 
A and B are ~ 0.25 and are smaller than in the other regions. These rati os coincide with 
those of H II regions in nearby disk galaxies (IMoustakas fc Kennicuttl 120061 ) . Therefore, we 
conclude that these reg ions are H II regions, w hich is why we call them H II regions A and 
B. As for H II region B, McCarthy et al.l (119871 ) observed a part of this region and suggested 
that it is an H II region. Our data show fairly uniform distribution of the [N II] /Ha ratio 
of ~ 0.25 over the whole H II region B including its three clumps. The [N Il]/Ha ratio of 
filament C is similar to those of H II regions A and B, rather than those of the galactic wind 
or starburst regions A and B. Therefore, this region is a star-forming region located on the 
same line of sight as an edge of the galactic wind although the region looks as if it is a part 
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qf the galactic wind. In Figure [2] (i) we compare our Ha map with the CO intensity map 
( IPaglione et al.ll2004l ). We find that the po sition of filament C coincides with that o f a strong 



CO emission elongated along the x<i orbit (IBinnev et al. 



1991 



Athanassoula Il992 ), which is 



made by and is oriented perpendicular to the bar (jPeng et al. 1996 ; Paglione et al. 2004 ). 



Hence, star formation in filament C occurs due to this £2 orbit. The [S II]AA6716,6731/Ha 
ratio map is globally similar to the [N II] /Ha ratio map: e.g., the [S II] /Ha ratios of filaments 
A and B are the largest among regions while those of H II regions A and B are the smallest. 
The [S II] /Ha ratios remain large outside of the filaments while the [N II] /Ha ratios become 
smaller. This may be because the metallicity of the galactic winds is larger than that outside 
of the filaments (see Section H~2l) . 

Figure 1(f) shows the [N II]A6583/[S II]AA6716,6731 ratio map. The [N II] /[S II] ratios 
in the central region and the galactic wind are almost uniform and are ~ 1.5, while the Ha 
fluxes vary. The ratio of filament C is small and is ~ 1.0. Two large [N II] /[S II] regions, 
whose ratios are greater than 2.0, are found to the north of starburst region A (Figure H]). 
Their positions coincide with those of star clusters found in HST observation (the bright 
blob and spot i: IWatson et al.lll996l ). The cause of the large [N II] /[S II] ratio is probably 
the nitrogen enrichment due to star clusters. We will discuss this in detail in Section 14.21 

Figure [2] (g) shows the Ha equivalent width map calculated from the Ha and continuum 
images. The equivalent width of H II region A is 95 A on average and larger than 150 A at 
its center, which is the largest among regions. Starburst regions A and B have local maxima 
of the equivalent width in the central region: 88 A and 64 A, respectively. Three clumps of 
H II region B have about the same equivalent widths with each other, ranging from 36 A to 
39 A. When the underlying diffuse continuum emission is subtracted, the equivalent widths 
of H II regions A and B are ~ 300-500 A (Table [1]). As diffuse continuum emission regions, 
we selected regions that are near the H II regions and do not include extra emission. This 
estimation of the Ha equivalent widths has an uncertainty of a factor 2 at most, mainly 
caused by the uncertainty of determination of this diffuse continuum level. 

The aperture, the Ha flux, the line ratios, and the equivalent width of each region are 
listed in Table [TJ The star formation rate is estimated from the Ha flux and also is listed 
in the table. As an example, we adopted the 6".0 diameter circular aperture in measuring 
values of H II region A. The measured Ha flux of H II region A is 3.5 xl0~ 13 erg cm -2 
s -1 , and its Ha luminosity is 5.1xl0 38 erg s _1 . From this Ha luminosity, the number of 
ionizing photons is esti mated to be 3.7x l0 50 s _1 and the star formation rate is estimated to 



be 4.0xl0~ 3 M yr" 1 flKennicuttlll998l ). 
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DISCUSSION 



4.1. Filament B 



At the tip of filament B, the [N Il]/Ha ratio reaches about 1.5 (see Section r3.2l) . An 
ordinary photoionization process in H II regions cannot produce thi s large [N II] /Ha ratio o f 
filament B, while shock ionization/excitation is plausible. Although ISchulz fe Wegnerl (Il992l ) 
produced the [N II] /Ha ratio larger than 1.0 with photoionization model, th is model used 



Kewlev et al. 



Kewley et al. 



a sing le stellar temperature for ionizing stars. The photoionization model by 
(2001) is more reliable, because they took into consideration stellar population. 
( 120011 ) placed a theoretical upper limit for starburst models on the optical diagnostic dia- 
grams (their equation (5) and Figure 16(a)), which show that it is difficult to exceed 1.0 in 
the [N II] /Ha ratio with starbursts even though they took into consideration the effects of 
the nitrogen enrichment. No combination of parameters, such as star formation history and 
ionization param eter, can generate a the oretical point above this fold. This limit is rather 
conservative, and lKauffmann et all (120031 ) placed a demarcation line at lower [N II] /Ha ratio 
from observational data. Photoionization may somewhat contribute near the central region, 
where the [N II] /Ha ratio is smaller than 1.0. 



According to the shock model calculated by lDopita fe Sutherland! (119951 ). the [N II] /Ha 
ratio reaches to 1.55, in the case of the solar-abundance interstellar gas, the magnetic pa- 
rameter B/n 1 ! 2 of 2 fiG cm 3 / 2 , which is the intermediate value among their models, and 
shock velocity of 500 km s _1 . This model can explain the largest [N II] /Ha ratio in our data, 
~ 1.5. By using the [S II]AA6716,6731/Ha ratio map and the [N II]/Ha ratio map, it is 
possible to discuss the possibility of nitrogen enrichment. Figur e [5] shows the [N II] /Ha rati o 
vs the [S II] /Ha ratio diagram of a shock model (shock-only: iDopita fe Sutherland! Il995l ). 
The observed ratios (squares and triangles) plotted above the model curves indicate that the 
N/O ratio is about 1.8 times the solar value when we assume that the magnetic parameter 
is equal to 2 and there is no contribution of photoionization for these ratios. We used the 
Anders fe Grevessel (1989) abundance set, because photoioniza tion and shock models that 
we r eferred to (IKewlev et all l200ll ; IDopita fe Sutherland! Il995l ) use it. If new abundance 
set (IGrevesse et all 120071 ) is used and if the model flux is proportional to the amount of 
each element, the N/O ratio will be 2.5 times the solar value. Both the [N II]/Ha and the 
[S II] /Ha ratios are larger in further positions from the galactic disk in filament B. This 
means that the shock velocity or the degree of contribution of photoionization varies among 
these positions. 

We compare the [N II] /Ha ratio of the galactic wind in NGC 253 with those in other 
galaxies. The [N II] /Ha ratio of the galactic wind in M82, which is one of the most famous 
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edge-on galaxies, is mainly less tha n 1.0, and photoionization by starb urst is considered to 
be an ionization/excitation source (jShopbell &: Bland-Hawthorn! Il998l ). There may be an 
influence of shock since the [N II] /Ha ratio tends to increase far from the galactic disk. In 
NGC 1482, in cont rast, the [N III /Ha ratio is larger than 1.0 in the bulk of the galactic wind 
and reaches to 2.0 (jVeilleux fc Rupkdl2002l ). This indicates shock as the ionization/excitation 
source. The [N II] /Ha ratio of galactic wind in NGC 253 is between the ratios observed 
for galactic winds in M82 and NG C 1482. The typical [N II I /Ha ratios are 0.3-0.6 in M82 
JShopbell & Bland-Hawtho7nlll998h and 1.0-2.3 in NGC 1482 (jveilleux fc Rupke!l2002h . while 
it is 0.8-1.5 in NGC 253. 

The shape of galactic wind and the line rat io maps are non-axisymmetric abo ut the 
galactic minor axis, which is also seen in M82 (jShopbell fc Bland-Hawthorn! Il998l ). The 
non-axisymmetric structure in NGC 253 is also fou nd in the H2 emiss ion, which suggests the 
inhomogeneous nature of the interstellar medium (ISugai et al.l 120031 ). Our [N II] /Ha ratio 
map shows that filament B has a larger ratio than filament A. This may occur if filament A 
is more affected by photoionization compared with filament B, or if the mass or the velocity 
of the wind gas varies from region to region. 



4.2. Nitrogen enrichment 



Starburst regions A and B emit the Ha and the continuum strongly and have low 
[N II] /Ha rat ios. Their [N II] /Ha ratios a re larger than those of H II regions in nearby 
disk galaxies ( jMoustakas fc Kennicuttll2006l ). and are about twice as large as those of H II 
regions A and B in NGC 253 itself (Figure [2] (d)). These large [N II] /Ha ratios are probably 
caused by nitrogen enrichment, as discussed in the following. Nitrog en enrichment rela tive to 
oxygen in metal rich environment was found in a statistical work by lLiang et al.l (120061 ). The 
nitrogen enrichment can easily explain the observed large [N II] /Ha as well as [N II] /[S II] 
ratio although distorted ionizing continuum, which mimicks the situation where neutra l 
matter lies between the ionizing stars and the emission-line clouds ( jSchulz &: Wegnerlll992l ). 
or shock may also increase the [N II] /Ha ratio. 



Nitrogen enrichment occurs in the following process (e.g.. |Pagellll997l ; ISugai et al.ll2004l ). 
There are two kinds of elements produced in stars, i.e., the primary element and the sec- 
ondary element. The primary element, such as oxygen and sulfur, is produced from hydrogen 
and helium. The ratio of the produced amount of each element (e.g., sulfur to oxygen) is 
constant and insensitive to the metallicity at the star birth. On the other hand the sec- 
ondary element, such as nitrogen, is produced also in the CNO cycle, and is spread into the 
interstellar medium with the stellar wind. The produced amount increases with an increase 
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in the metallicity at the star birth. Therefore, in a metal rich environment, the ratio of 
the produced amount of nitrogen (secondary element) to that of sulfur (primary element) 
becomes larger, and the [N II] /[S II] ratio becomes large. The relationship between the 



[N II] / [S II] ratio and the metallicity of galaxies is found in SDSS observations (ILiang et al. 



2006|). 



In the [N II] /[S II] ratio map (Figure [2] (f)), the [N II] /[S II] ratios of starburst regions 
A and B are larger than those of H II regions A and B and filament C. This indicates that 
starburst regions A and B have higher metallicity. The [N II] /[S II] ratios in the bright blob 
and spot i (hereafter, we call it blob/spot region) are the largest (Figur elU), so that the degree 
of nit rogen enrichment is extremely high. The de-excitation of [S II] (lOsterbrock fc Ferland 



20061 ) may cause this large [N II] /[S II] ratio. The electron densities ar e estimated to be 
no mo re than ~ 900 cm -3 from [S II]A6716/[S II] A6731 ratios observed by lSchulz &: Wegner 
([1992]). This ratio is at the minimum (~ 0.9) around the blob/spot region in the whole 
field of view of their spectroscopic observations. These electron densities are small and the 
de-excitation reduces the [ S II] flux by no m ore than 16 %. The bright blob and spot i are 
bright in J, H, and K band (jSams et al. I1994T) . and are two out of four clusters defined in the 
optical continuum in HST observation (IWatson et al.lll996l ). The bright blob is the brightest 
cluster among the four. Spot i is the second brightest in the Ha flux. We estimate the 
metallicity of blob/spot region from the [N II] /[S II] rati o, 2.19. This ratio is at the largest 
end of the distribution for starburst galaxies analyzed by ILiang et al.l (120061 ). Assuming the 
rel ation between the [N II] /[S II] ratio and the metallicity in starburst galaxies (Figure 3(d) 
bv lLiang et~al~1l2006h . the metallicity (12+log(0/H)) of blob/spot region is der ived as 9.17 
± 0.08 (lcr). Using the relation between log(0/H) and log(N/0) (Figure 8 by ILiang et al. 



20061 ). the log(N/0) abundance of blob /spot region is estimate d to be -0.50 ± 0.13 (lcr). 



Here the scatter in N/O to O/H in Figure 8 bv lLiang et al.l (12006!) was estim ated to be about 
0.1 dex (lcr) by eye-inspection. When we use the lAnders fc Grevessd (Il989l ) abundance set, 
whose log(N/0) ratio is -0.88, the N/O ratio of blob/spot region is about 2.4 times larger 



than t he solar value. This value of 2.4 does not change even if we use the iGrevesse et al. 
(120071 ) abundance set, because its N/O ratio is equal to that of the lAnders &: Grevesse Jl989h 
abundance set. This value is similar to that in the galactic wind (Section 14. II) . although it is 
difficult to estimate the abundance precisely. To explain this similarity, it is plausible that 
the gas of the galactic wind is chemically evolved in starburst regions and then is blown out. 
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4.3. Extent of the galactic wind 



In H II region B, three clumps are found in our Ha map. The northern and western 
clumps are on the spiral arm, while the southeastern one is not. Soft X-ray (0.3-2.0 keV) 
is detected around H II region B, including the three clumps (Figures [2] (h) and [3] (b)). 
Therefore, it is difficult to determine from only the Ha map whether the southeastern clump 
is a part of the galactic wind or an H II region on the galactic disk with an X-ray emitting 
gas accidentally located on the same line of sight. This issue is important when considering 
the structure of the galactic wind, i.e., whether or not the tip of the galactic wind is cooled 
down enough to emit the Ha. Although H II region B was observe d and considered to be 
H II regions through a long-slit spectroscopy ([McCarthy et al.l 119871 ). the line ratio of each 
clump was not discussed. In our data, the southeastern clump, as well as the northern and 
the weste rn ones, have the same [N II] /H a and [S II] /Ha ratios as H II regions in nearby disk 
galaxies (IMoustakas fc Kennicuttl 120061 ). Moreover, the equivalent widths of three clumps 
are similar: ~ 40 A (Figure [2] (g)). Consequently, we consider that the three clumps of H II 
region B are not a part of the galactic wind but are H II regions on the galactic disk. 



4.4. Kinetic energy of the galactic wind 



We derive the properties of the galactic wind, given that the region where the [N II] /Ha 
ratio is greater than 1.0 (or the [S Il]/Ha ratio is greater than 0.5) is the galactic wind 
dVeilleux etaDbood : iBland-Hawthorn et al.l 120071 ). First we estimate the Ha luminosity of 
the wind. When the [N II] /Ha criterion is used, the Ha flux of the galactic wind is equal 
to 1.6 xlO -13 erg cm -2 s -1 , and its Ha luminosity is 2.3 xlO 38 erg s _1 when 3.5 Mpc is 



used as a distance to NGC 253 (IRekola et al.ll2005l ). When the [S II] /Ha criterion is used 



instead, the Ha flux is equal to 8.9xl0~ 13 erg cm~ 2 s _1 and the Ha luminosity is 1.3 xlO 39 
erg s -1 . Hereafter, we adopt the [S II] /Ha criterion, because it is less affected by nitrogen 
enrichment (Section 14.11 14.21) compared with the [N II] /Ha criterion. Bec ause the total 



Ho+[N II]AA6548,6583 luminosity of NGC 253 is equal to 1.1 xlO 41 erg s" 1 flHoopes et al.l 



19961 ). the Ha luminosity of the galactic wind is only from 0.5% to 5% of the total luminosity. 
The actual fraction of the luminosity of the galactic wind may be even smaller, because 
extinction is not corre cted in the above estimate although the central region of NGC 253 
is heavily dusty (e.g., ISams et al.l ll994T ) . Secondly, we estimate the mass of ionized gas 
in the galactic wind. From th e [S II]A6716/[S Il]A6731 ratio of 0.9 to 1.2 in the central 
region JSchulz fc Wegnerlll992h. t he ele ctron density is estimated to be from 100 cm to 



1000 cm" ~ 3 JOsterbrock & Ferlandlhood ). Using the electron density of 100n e 2 cm 3 , where 



n e 2 is normalized to 100 cm 3 , and the recombination rate of Ha of 8.7x10 



14 S —1 

cm s 
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( IQsterbrock fc Ferlandll2006l ). the Ha emitting gas mass of the galactic wind is estimated to 



be 8.6xl0 4 n~2 M Q . We assumed that the northwest side of the galactic wind, which is not 
detected in our observations due to heavy extinction, actually exists as well as the southeast 
side of it (Section 13. II) . For simplicity, we assumed the symmetric shape of the galactic wind 
in order to roughly estimate the volume filling factor, although it is non-axisymmetric about 
the galactic minor axis. The assumed shape of the galactic wind is a truncated cone structure 
whose opening angle is 30° and diameter of the base and slant height are 500 pc and 700 
pc, respectively. These values are measured through the Ha map. Then, the volume filling 



factor of the galactic wind is estimated to be 4.2 x 10" 



n 



e,2- 



The limb-brightening in the Ha 



and the X-ray (Section 13. ip and the small filling factor suggest tha t the galactic wind has a 
hollow structure. Finally, using the outflow velocity of 390 km s -1 ( ISchulz &: Wegnerlll992l ). 



the kinetic energy of the galactic wind is estimated to be 1.3xl0 53 n e \ erg. The dynamical 
time scale Td yn derived from the length of the filament B (~ 700pc) and the outflow velocity 



(390 km b _1 ; 
average. 



is ~ 1.8 Myr, and the kinetic energy is increased by 2.3x10 



n 



e,2 



erg s 



on 



In order to consider whether the energy input from supernovae can supply the kinetic 
energy of the galactic wind, we estimate first the Ha luminosity of the bright blob, which is 
the brightest region in the central region of NGC 253. Because of heavy extinction in the 
central region of NGC 2 53, we have correct ed the effect of extinction using the ratio of the 
Ha flux to t he Br7 flux (jForbes et al.lll993T ). We assumed the intrinsic Br7/Ha flux ratio of 



9.80xl0~ 3 flOsterbrock fc Ferland 



A(2.2/im)=0.15A y flWhitford 



1958 



200 



^) and the extinction law of A(6563A)=0.79Ay and 



Miller fc Mathewslll972j^T he Ha flux of the bright blob 



in 2"aperture, which is the same as used by lForbes et al.l (119931 ). is 5.5x10 14 erg cm 2 s 1 , 



and Ay is estimated to 6.9 mag. The extinction-corrected Ha flux is 8.3x10 12 erg cm 2 s 1 , 
and the Ha lu minosity is 1.2 xl O 40 erg s _1 , which corresponds to the star formation rate of 



0.096 M Q yr _1 (jKennicuttlll998l ) . When we assume that the star formation rate of the bright 
blob is constant, the supplied kinetic energy through supernovae and stellar winds from thi s 
region in 1.8 Myr (= Tdyn) is estimated to be 5.6 xlO 53 erg (S tarburst99: lLeitherer et al.lll999l ). 
We assumed Salpeter initial mass function (jSalpeterlll955l ). M up =lOOM and Mi ow =lM Q as 
cut-off masses, metallicity Z=0.020, and that the galactic wind be gan to blow when th e 
first supernova occurred in the bright blob (3.6 Myr: Figure 114 in iLeitherer et al.lll999l ). 
Therefore, the energy input from the bright blob is much greater than the kinetic energy of 
the galactic wind. This suggests that starbursts are the main energy source of the galactic 
wind in NGC 253. 



W e compare the properties of the ga lactic winds in M82 (IShopbell fc Bland-Hawthorn 
1998T) . NGC 1482 JVeilleux fc Rupke!l2002h . and NGC 253, although the de finition of galactic 
wind regions in M82 is quite different from that in NGC 253 and NGC 1482. IShopbell fe Bland-Hawthorn 
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( 119981 ) defined the galactic wind regions by simply excluding with in approximately 8"of the 
disk (+ halo component), while we and IVeilleux fc Rupkel (120021 ) used the criterion of the 
[S II] /Ha and [N II] /Ha ratio, respectively. The kinetic energy of the galactic wind in NGC 
253 (1.3 xlO 53 erg) is found to be smaller than those in M82 (2.1xl0 55 erg) and in NGC 1482 
(2xl0 53 erg). There are three possibilities for the difference in the energies of the galactic 
winds: the masses, the velocities, and the dynamical timescales. The mass in NGC 1482 
is about 4.2 times, and the mass in M82 is about 67 times as large as that in NGC 253 
(8.6xlO 4 M ). The outflow velocities and the dynamical timescales do not vary drastically 
among these galaxies: 390 km s _1 and 1.8 Myr in NGC 253, 250 km s -1 and 6 Myr in NGC 
1482, and 600 km s _1 and 3.4 Myr in M82, respectively. Therefore, the difference of the 
outflow mass in each galaxy, which may depend on the strength of starburst in the central 
region, probably dominantly makes the variety of the kinetic energy of the galactic wind, 
although the combination of the differences of the mass, the velocity, and the dynamical 
time scale will determine the exact value. 



5. SUMMARY 

We have observed the central region of NGC 253 with the Fabry-Perot mode of the 
Kyoto3DII mounted on the Subaru Telescope. We have presented the continuum and the Ha 
images, and the [N II]A6583/Ha, [S II]AA6716,6731/Ha, and [N II]A6583/[S II]AA6716,6731 
ratio maps. The [N II]/Ha of filament B is larger than 1.0, which indicates that this region is 
ionized/excited by shock. Non-axisymmetry in the galactic wind, which is also seen in M82, 
is found. The shapes of the galactic wind and the line ratio maps are non-axisymmetric about 
the galactic minor axis: e.g., the [N II] /Ha of filament A is smaller than that of filament 
B. The [N II] /Ha map enables us to isolate filament C from the galactic wind. Using the 
[S II] /Ha map, we extract the galactic wind region, and estimate the mass and the kinetic 
energy of the galactic wind: 8.6xl0 4 M Q and 1.3 xlO 53 erg, respectively. The mass and the 
kinetic energy of the galactic wind in NGC 253 are smaller than those in M82 and NGC 
1482. The filling factor of the galactic wind is much smaller than unity, which is consistent 
with that the galactic wind has a hollow structure. The energy input from supernovae in 
the bright blob can sufficiently supply the energy of the galactic wind. From the [N II] /[S II] 
ratio map, the interstellar gas in the central region and in the galactic wind is found to have 
a high nitrogen enrichment. The bright blob and spot i, which are stellar clusters found in 
HST observations, have especially high [N II] /[S II] ratios due to heavy nitrogen enrichment. 
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Fig. 1. — Good (left panel) and poor examples (right panel) of fitting of the Ho;+[N II] set. 
The flux is shown as a function of wavelength in units of A. The observed values are plotted 
in filled circles. The dashed line and dotted line represent the fluxes from the Ha and the 
[N II]AA6548,6583, respectively, as a function of the central wavelength of the transmission 
curve of etalon. The solid line is sum of the dashed line and dotted line. 
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Fig. 2. — The continuum and Ha image, and the line ratio maps of the central region of 
NGC 253. North is up and east is to the left. The green cross at t he center of each figure 



repre sents the center of NGC 253 derived from the radio continuum (jUlvestad fc Antonucci 



19971 ). The bars in all figures at lower left represent the length of 100 pc. (a) The continuum 



1 arcsec 2 A 1 . (b) The deep and (c) shallow Ha image 



image in units of 10 -16 erg cm -2 s" 
in units of 10 -16 erg cm -2 s -1 arcsec -2 . (d) The [N II] /Ha line ratio map. The background 
noise level in the Ha flux (la) is 1.9 x 10 -16 erg cm -2 s -1 arcsec" 1 . The green contours are 
4cr, 8cr, 16a, 32<t, 64<t, and 128cr of the Ha flux. The region where the Ha flux is weaker than 
3a (= 5.7 x 10 -16 erg cm -2 s -1 arcsec -2 ) is shown in black. The blue squares represent the 
regions where the [N II] /Ha and [S II] /Ha ratios are measured and are plotted with green 
squares in FigureEl (e) The [S II]/Ha line ratio map. (f) The [N II] /[S II] line ratio map. 
(g) The Ha equivalent width map in units of A. (h) The continuum subtracted Ha image in 
units of 10 -16 erg cm -2 s -1 arcsec -2 . The contours show the intensity of soft X-ray (0.3-2.0 
keV: IStrickland et al.ll2000l ). (i) The continuum subtracted Ha imag e in units of 10 -16 er g 
cm -2 s -1 arcsec -2 . The contours show the intensity of CO J = 1 — > (IPaglione et al.ll2004f ). 
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3. — (a) The continuum subtracted Ho images of NGC 253 in units of erg cm 2 s 1 
2 . North is up and east is to the left. The green cross at t he center of this figure 



Fig 

arcsec 

repre sents the center of NGC 253 derived from the radio continuum (jUlvestad &: Antonucci 



19971 ). The bar at lower left represents the length of 1 kpc. The region where the Ha flux is 
weaker than 3a (= 5.7 x 10 -16 erg cm -2 s -1 arcsec -2 ) is shown in black. The bright point 
at the lower right edge is a ghost of the central region, (b) The continuum subtracted Ha 



images of NGC 253 in units of erg cm 2 s 



_1 arcsec 2 . 



soft X-ray (0.3-2.0 keV: IStrickland et alJl2000h . 



The contours show the intensity of 
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Fig. 4. — The [N II]A6583/[S II]AA6716,6731 ratio map in the central region. North is up 
and east is to the left. The cross at th e center of this figure repres ents the center of NGC 
253 derived from the radio continuum (jUfvestad Antonuccil 119971 ). The bar at lower left 
represents the length of 100 pc. Four plus signs represent the positions of star clusters found 
by HST, which are spot a, spot n, spot i, and the bright blob starting from the left. The 
region where the Ha flux is weaker than 3a (= 5.7 x 10~ 16 erg cm -2 s _1 arcsec~ 2 ) is shown 
in black. [See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 5. — The [N II] /Ha ratio y s the [S II] /Ha ratio diagram of the shock model (shock- 
only: iDopita fc Sutherlandlll995l ). The numbers below the model curves represent the shock 
velocities in units of km s , and the numbers on the curves represent the magnetic pa- 
rameters (B/n l l 2 ) in units of /iG cm 3 '' 2 . The effect of a factor of 2 increase in the N/O 
ratio is shown by the arrow. Filled circles are for LINERs while crosses are for radio galaxies 
JPopita fc SutherlandlE995l l. 7 our squares represent the line ratios observed at filament B in 
NGC 253 (see Figure [2] (d) for their positions) while five triangles represent those in typical 
regions of the galactic wind. The square drawn furthest left shows the line ratios closest to 
the galactic disk, with each subsequent square to the right lying at increasing distance. [See 
the electronic edition of the Journal for a color version of this figure] 



Table 1. Fluxes, star formation rates, and line ratios of each region 



Region 


Aperture 


Flux(Ha) a 


SFR b 


EW(Ha) c 


[N II]/Ha d 


[S II]/Ha c 


[N II] /[S II] d ' c 


starburst region A 


3". 5 


1.9 


2.2 


88 


0.54 


0.28 


1.93 


starburst region B 


4".0 


2.5 


2.9 


64 


0.51 


0.35 


1.45 


filament C 


11". 54 x 3" .47 


3.0 


3.5 


42 


0.38 


0.32 


1.16 


H II region A 


6".0 


3.5 


4.0 


95(482) 


0.25 


0.17 


1.45 


H II region B(N) 


4".0 


0.34 


0.39 


36(316) 


0.24 


0.22 


1.06 


H II region B(W) 


6".0 


0.72 


0.83 


36(371) 


0.26 


0.24 


1.09 


H II region B(SE) 


8".0 


1.1 


1.3 


39(282) 


0.26 


0.26 


0.98 



Note. — The errors for flux(Ha), EW(Ha), [N II] /Ha, [S II]/Ha, and [N II] /[S II] are less than 30%, a factor 
of 2, 30%, 10%, and 30%, respectively. 

a The unit of flux is 10~ 13 erg cm -2 s _1 . 



b The star formation rates are estimated from the Ha fluxes using the relation obtained by iKennicuttl (119981 ) . 
The unit of star formation rate is 10 -3 M yr _1 . 

c The unit of equivalent width is A. The value shown in parentheses is the equivalent width derived by 
subtracting the underlying diffuse continuum from the continuum emission. 

d [N II] means the [N II]A6583 line. 

e [S II] means the [S II]AA6716,6731 lines. 



